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For the observation of Bose-Einstein condensation, ex- 
citons in cuprous oxide are regarded as promising can- 
didates due to their large binding energy and long life- 
time. High particle densities may be achieved by entrap- 
ment in a stress induced potential. We consider a multi- 
component gas of interacting para- and orthoexcitons in 



cuprous oxide confined in a three-dimensional potential 
trap. Based on the Hartree-Fock-Bogoliubov theory, we 
calculate density profiles as well as decay luminescence 
spectra which exhibit signatures of the separation of the 
Bose-condensed phases. 
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1 Introduction The theoretical framework for trapped 
dilute interacting bosonic gases is well known from the 
theory of atomic condensates [1 ,2,3,4]. First applications 
to excitonic systems exist, as well (5). Recent investi- 
gations in the framework of a mean-field formalism in 
local density approximation suggest distinct signatures 
of a condensate in the decay luminescence spectrum of 
the thermal excitons JS ,7 ,8|. Works on two-component 
systems 19 10i lTTlfT2l have shown that the occurrence of 
phase separation is closely tied to the proportions of inter- 
and intra-species interaction strengths. Yet, the description 
of exciton-exciton interaction is a long-standing problem. 
Experimental results and theoretical predictions for the in- 
teraction strengths vary within an order of magnitude lfl3l 
I5l fl4ll . As an example, we present numerical results for the 
densities and the spatially resolved luminescence spectra 
of the three component system of excitons in cuprous ox- 
ide (Cu 2 0), i.e., para-, ortho(+)-, and orfho(— )excitons, 
trapped in a strain induced potential |15|. We show how 
spectral features may reveal phase separation, thereby 
yielding a minimum estimate of the relative strength of the 
mutual interactions. 

2 Multicomponent exciton systems 

2.1 Thermodynamics We consider a X-component 
exciton gas in second quantization, starting from the 



Hamiltonian in the grand canonical ensemble: 

i K r 

+ 2 ^ / d 8 '^.^'.*)^^.*)^^*)^-^,*) ,(1) 

where Vi represents the external (trap) potentials and //$ 
the chemical potentials. We assume a contact potential for 
the exciton-exciton interaction, with the intra- and inter- 
species interaction strengths hij — 2irh 2 (M^ 1 +M~ 1 )afj 
given by the respective s-wave scattering lengths afj . 

The Bose field operators ^ (r , t) are decomposed in the 
usual fashion, ipi(r, t) = <£j(r) + V>i(r, t), with the conden- 
sate wave functions <Pi(r) = (^(r,i)) = (ipi(r)) and the 
operators of the thermal excitons if>i(r, t). The Heisenberg 
equations of motion \hdtipi — [ipi,W] result in 2K cou- 
pled equations (arguments dropped for brevity): the Gross- 
Pitaevskii equations (GPE) for the condensates, 

0= y- +Vi- ii, + ha (n lt + n u ) /',,/',, J '/', 

+ hatha®* + h l0 (n,,<I>, + ///„'/':) , (2) 
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and the equations of motion for the thermal excitons, 



2Mi 



+ Vi- /ii + Zhan^ +^ h t 



Hi 



(3) 



Here = $*<l>i + n.y, = + my, with the aver- 
ages riij = (ipjipj) and my = (ipiipj). For simplicity, we 
neglect all non-diagonal averages, i.e., the last terms on the 
r.h.s of Q and ([3]) and subsequently obtain effective one- 
component equations with mean field contributions from 
the respective other components. Because the extension 
of the potential trap is large compared to the thermal de- 
Broglie wavelength of the excitons, we apply a local den- 
sity approximation to (^J, setting V 2 — > — |k| 2 with a wave 
vector k. For the same reason, we apply the Thomas-Fermi 
approximation to the GPE, thus neglecting the kinetic en- 
ergy term in pi. With the above simplifications, Eq. Q is 
solved by a Bogoliubov transformation 



ipi = 



£ [ Ml (kK(k)e- i ^ k >'/' 1 + «?(k)al(k)e 



iEi(k)t/h 



(4) 

The densities nj = nu of thermally excited excitons are 
given by 



d 3 k 

8tt3 



L,(k. r) / . .. 1 

4^(^(k,r)) + 2 



xete(k,r) 



(5) 



with ub{E) = [exp(i?/fcsT) — 1]~ . To guarantee gap- 
less excitation spectra Ei, we neglect all anomalous aver- 
ages fhu (Popov approximation) and obtain 



£ i (k,r) = v /L i (k,r)2-(/ lii n?(r)) 2 
?i 2 £; 2 



(6) 



Li(k,r): 



2M, 



Vi(r) - ^ + 2h ll n l (r) + h tj nj(r) , 

(7) 



with = |^| 2 and = nu — nf + nf. 

From the simplified GPEs, the condensate densities fol- 
low as 



<(r) 



(8) 



if this expression is non-negative, and nf(r) = oth- 
erwise. Equations |5]l to ([8]l have to be solved self- 
consistently. Although they look similar to the one- 
component case, the coupling between the components 
appears in and n\. 

2.2 Luminescence spectrum Excitons decay by 
emitting photons. We apply a local approximation to the 



emission spectrum, which is determined by the excitonic 
spectral function A(k, ui) B16II17I : 

Ii( r ,u) oc 27r|5<(k = 0)| 2 <5(W - ^K(r) (9) 
+ ^2\Si(k)\ 2 n B (hL)' - iM)Ai(r,\s, hJ - m) , 



with the exciton-photon coupling S'i(k). The spectral func- 
tion is given by the Bogoliubov amplitudes Ui and Uj, and 
by the quasiparticle spectrum in dSJ: 



Ai(r, k, w) = 2nh 



(10) 



W 2 (k,r)<5(^ + ^(k,r)) 



In CU2O, the decay of orthoexcitons takes place via mo- 
mentum supplying phonons, such that all exciton states k 
participate and uj' = us — E g x /K — i^phonon , with E g x be- 
ing the excitonic band gap. The paraexcitons decay with- 
out phonons, i.e., ui' = ui — E g x/h. Here, energy and 
momentum conservation only allow for processes, where 
the wave vectors of excitons and photons are equal, so that 
S(k) = S S(k - k ), with |k | = E gX n/hc. Therefore, 
the condensate does not contribute to the zero-phonon pro- 
cess. Despite this, there are indirect signatures of the con- 
densate in the spatially resolved luminescence spectrum of 
the non-condensed excitons |6] . 

In a typical experiment, one images a small stripe of 
width 2Ax elongated along the z-direction onto the en- 
trance slit of a spectrograph. Thus, by integrating ([9]) over 
the x- and y-directions perpendicular to z, we obtain the 
spatially resolved spectrum. Moreover, we account for the 
finite spectral resolution A by convoluting the spectral in- 
tensity with a slit function of the shape exp{— (ui/A) 2 }. 
For our calculations, we use values of A — 25 \xeV for 
the spectral resolution and Ax = 25 (i.m for the entrance 
slit of the spectrograph being typical for a triple high- 
resolution spectrograph used in the current experiments 
which are underway (TSJ. 

3 Results We calculated the density distributions 
and the luminescence spectra of the trapped excitons for 
three different choices of the interaction strength between 

ortho(+)- and ortho(— )excitons (h^ ). The interaction 

strengths are calculated from the scattering lengths given 

in lfT3l (h++ = h = 0.71h pp , h p+ = h p - = 0.33h pp , 

and h-\ = 1.77h pp with h pp — 5.4 x 1CP 4 |xeV(i.m 3 ). 

Taking the value for h-\ given by Shumway and Ceper- 

ley llT3ll as hsc we chose: (i) h-\ = hsc/5; (ii) 

= /isc/3; and (hi) h + - = hsc while keeping 
the others fixed. The used trap potentials Vi are fitted to 
experimental data taken from [18|. The minimum of the 
paraexciton trap is —1981 |J.eV and the minimum of the 
orthoexciton trap is —8158 \ieV. Note that we neglect the 
difference in paraexciton and orthoexciton mass due to the 
k-dependent exchange interaction 1 19 1. 
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The results of the numerical calculations are shown in 
Fig.[T] The densities are visualized as contour plots includ- 
ing thermal and condensed excitons using cylindrical coor- 
dinates (/0,z). The densities are independent of tp due to the 
trap geometry. The thermal excitons form the darker rim 
of the bright condensate spots. The second rows show the 
luminescence spectrum of the paraexcitons and the com- 
bined spectrum of the orthoexcitons. The latter contains 
the condensate signals, i.e., the first term on the r.h.s of (|9j, 
assuming S(k) is a k-independent constant. 

At a temperature of T = 0.7 K, all three exciton 
species show Bose-Einstein condensation. The condition 
for phase separation of the orthoexcitons is approximately 
h\_ > h ++ h [9] and becomes in the considered sys- 
tem > h ++ since h ++ = h In case (i) where h-\ 

is considerably smaller than h ++ the ortho-condensates are 
nearly completely mixed. Due to the different minimum 
positions of the respective external potentials, the thermal 
paraexcitons are pushed aside by the combined orthoexci- 
ton densities. The spectrum of the orthoexcitons is domi- 
nated by two condensate peaks at the chemical potentials 
^_ = -8013 u.eV and /i+ = -7990 u.eV. In case of 
the paraexcitons, only the zero-phonon spectral line con- 
tributes. It reveals the existence of a condensate via the flat 
bottom at e - fi p = GEM). 

Increasing to hsc/3 yields the results shown in 

case (ii). Here is slightly smaller than h ++ and the 

condensates are still mixed in a wide area. However, a 
starting of the separation can be observed. Nevertheless, 
the spectra for the para- and orthoexcitons as well as the 
paraexciton density are not changed qualitatively with re- 
spect to case (i). 

In case (iii) the condition for phase separation, > 

h ++ , is fulfilled and the ortho-condensates form a ball- 
and-shell structure with finite overlap. Incidentally, the dif- 
ference in their chemical potentials is smaller than the 
spectral resolution so that the combined ortho-spectrum 
gives no evidence of the phase separation. However, in the 
region of overlapping ortho-condensates we find a notice- 
able depletion of paraexcitons, which results in a W-shaped 
distortion of the para-spectrum when compared to case (i) 
or (ii). 

When the interaction strength h-\ is further increased, 

no qualitative changes with respect to case (iii) are found. 
For T > there are no pure (+) or (— ) phases lfT2l and the 
spectral features described in case (iii) remain. 

4 Conclusion and outlook Our simulations for 
the experimentally relevant example cuprous oxide have 
shown that at finite temperatures a possible phase sepa- 
ration of excitonic condensates may not be reflected in 
their combined luminescence spectrum. However, in the 
case of cuprous oxide, the single spectral line of the di- 
rect paraexciton decay may consitute a rather sensitive 
probe for the spatial structure of the orthoexciton density 
distribution. Specific distortions in the paraexciton spec- 



trum would provide an experimental footprint of a phase 
separation of orthoexciton condensates. If detected, the in- 
teraction strength of the (+) and (— ) species would satisfy 

the relation h-\ > ft++. 

Omitting the Thomas-Fermi approximation or one 
of the other approximations used, might lead to results 
that quantitatively differ from the ones presented here. 
However, the above conclusions are mainly based on the 
existence of a finite overlap of separated condensates 
and should remain valid even if the spectra and densi- 
ties are calculated beyond the used approximations. Also 
the calculations for the luminescence spectrum could be 
enhanced using more advanced approaches 11201 . Further 
research is required to address these issues. 
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Figure 1 Three different subsets of exciton densities and luminescence spectra for para-, ortho(+)- and 
ortho(— )excitons for T = 0.7 K and a total particle number of N p = 10 10 paraexcitons and iV_ = N + = 10 10 
orthoexcitons in the trap. Each subset contains 5 images with the exciton densities in the first row (from left 
to right para-, ortho(+)- and ortho(— )excitons) and the spectra for the para- (left) and orthoexcitons (right) in 

the second row. Taking given by Shumway and Ceperley lfl"3l as hsc the subsets represent from top to 

bottom: (i) = hsc /5; (ii) = hsc /3; and (iii) h^ = hsc while taking all other interaction strengths 

as given in lfT3ll . 
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